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ABSTRACT 

We present high-precision relative radial velocities of the very low-mass star VB 10 that were ob- 
tained over a time span of 0.61 yr as part of an ongoing search for planets around stars at the end 
of the main sequence. The radial velocities were measured from high-resolution near-infrared spectra 
obtained using the CRIRES instrument on the VLT with an ammonia gas cell. The typical internal 
precision of the measurements is 10ms . These data do not exhibit significant variability and are 
essentially constant at a level consistent with the measurement uncertainties. Therefore, we do not 
detect the radial velocity variations of VB 10 expected due to the presence of an orbiting giant planet 
similar to that recently proposed by Pravdo and Shaklan based on apparent astrometric perturba- 
tions. In addition, we do not confirm the ~ 1 kms -1 radial velocity variability of the star tentatively 
detected by Zapatero Osorio and colleagues with lower precision measurements. Our measurements 
rule out planets with M p > 3 Mj up and the orbital period and inclination suggested by Pravdo and 
Shaklan at better than 5 a confidence. We conclude that the planet detection claimed by Pravdo and 
Shaklan is spurious on the basis of this result. Although the outcome of this work is a non-detection, it 
illustrates the potential of using ammonia cell radial velocities to detect planets around very low-mass 
stars. 

Subject headings: planetary systems — stars: individual (GJ752B) 



1. INTRODUCTION 

Very low-mass stars are an interesting sample of po- 
tential planet hosts. They are useful for probing the 
correlatio n between gas giant planet frequency and stel- 
lar mass (|Endl et al.ll2006t I Johnson et all 12007ft . are the 
most numerous stars in the Galaxy (jHenrv et al.l f2006). 
and exhibit a larger dynamical response to orbiting plan- 
ets a nd have closer-in ha bitable zones than higher-mass 
stars (jKasting et al.lfl993l ). However, the search for plan- 
ets around these stars is hindered by their extreme in- 
trinsic faintness at visible wavelengths, and few sensitive 
searches for planets around the m have been carried out. 

Recently, iPravdo fc Shaklan! (|2009l hereafter PS09) 
have claimed the first detection of a planet around 
a nearby, very low-mass star based on apparent peri- 
odic astrometric perturbations. The star is VB 10 and 
the planet candidate is a giant planet with a mass 
M p = 6.4Mj U p and orbital period P = 0.744yr. This 
could also be the first astrometric discovery of an ex- 
oplanct. 

We have been observing VB 10 as part of an ongo- 
ing search for planets around the lowest-mass stars (i.e. 
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]VL<O.2M ) using high-preci sion near-infrared (NIR) 
radial velocity measurements (jBean et~aT1[2009l) . The 
advantages of the NIR for radial velocity measurements 
are that it is possible to obtain high-resolution and high- 
signal-to-noise spectra of very low-mass stars in these 
wavelengths with current instruments, and also that the 
affect of activity induc ed "jitter" is poten tially reduced 
relative to the visible (jReiners et al.1 120091 ) . We present 
here our radial velocity measurements of VB 10 obtained 
so far, and discuss them in the li g ht of PS09's claimed 
detection. IZapatero Osorio et all (|2009l hereafter Z09) 
have presented NIR radial velocities with a typical in- 
ternal precision of 300 ms" 1 that were obtained using 
NIRSPEC on Keck ([McLean et alJl998t l. These data ex- 
hibit variability on the level of ~ lkms , which would 
be consistent with the existence of a giant planet. Our 
radial velocity measurements, which are 30 times more 
precise than those in the Z09 study, also allow us to make 
a direct comparison to their results. 

2. RADIAL VELOCITY MEASUREMENTS 

We obtained spectroscopic observations of VB 10 for 
radial velocity measurement s using the high- resolution 
NIR spectrograph CRIRES (|Kaufl et al.ll2004l) . which is 
fed by the UT1 telescope of the Very Large Telescope 
(VLT) facility. As with all the measurements obtained 
as part of our CRIRES planet search, a glass cell filled 
with ammonia gas was placed in front of the spectrograph 
entrance to enable high-precision calibration of the in- 
strument response during analysis of the data. We used 
a 0'.'4 slit for the VB 10 observations, which corresponds 
to a nominal resolving power R = A/AA « 50000, and 
we recorded spectra in the CO bandhead region of the K- 
band (~ 2.3 /im). This is a region in which late-type stars 
exhibit useful lines for radial velocity measurements and 
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TABLE 1 
Radial Velocities Of VB 10 
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the ammonia cell also exhibits useful calibration lines. A 
total of 36.4 nm of the 48.2 nm of recorded wavelength 
coverage was utilized for the radial velocity measure- 
ments. 

Observations were obtained during 12 nights in four 
classical style runs between March and November 2009. 
Each visit consisted of a sequence (typically four) of 300 s 
exposures. We reduced the data and measured relative 
radial velocities from the extracted spectr a using the 
meth ods described by us in a previous paper (|Bean et al.l 
2009). Our radial velocity measurement algorithm is an 
ad aptation of the cano nical gas cell technique described 
bv lButler et al.l (|1996l ). and includes simultaneous mod- 
eling of the wavelength scale and instrument point spread 
function in each spectrum using the ammonia cell lines 
as a fiducial. We have previously demonstrated that our 
method can deliver long-term precisions of 5ms _1 for 
late M dwarfs when applied to data of sufficient qual- 
ity (i.e. signal-to-noise and resolution) obtained with 
CRIRES using the ammonia cell. We concentrated only 
on measuring and analyzing relative radial velocities for 
this study because our gas cell method enables us to 
reach velocity precisions on an internal scale that are 
roughly an order of magnitude better than we could ob- 
tain for velocities on an absolute scale. The higher preci- 
sion we obtain for relative as opposed to absolute veloci- 
ties yields the best sensitivity for the detection of unseen 
companions orbiting VB 10. 

The signal-to-noise ratios (SNR) of the reduced and 
extracted Id spectra of VB 10 are ~100pixel _1 . The 
expected photon-limited radial velocity precisions from 
these data are typically 13 ms -1 including calibration 
uncertainty, and the median of our estimated uncer- 
tainties on the radial velocity measurements is 22 ms -1 . 
Our estimated uncertainties are normally larger than 
the corresponding photon-limited precisions because our 
error estimation method also accounts for the imper- 
fect spectrum modeling during the measurement process. 
Previous observations of radial velocity standard stars 
have shown that our adopted measurement uncertain- 
ties are reasonable (down to the long-term noise floor of 
~5ms _1 ), and we have found that the determined un- 
certainties for our method are usually a fac tor of 1.5 - 
2.5 la rger than the photon-limited precisions (|Bean et al.l 
2009L That our estimated uncertainties on the velocities 
in this case are typically 1.7 times the expected photon- 
limited error is completely consistent with this previous 
result. Ultimately, our VB 10 radial velocities are pri- 
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Fig. 1. — Radial velocities for VB 10 measured using CRIRES 
at the VLT with an ammonia gas cell. Multiple exposures taken 
within up to 2.5 h have been binned to give each point. 

marily limited by the obtained SNR in both the spectra 
that were analyzed to give the measurements, and also 
the template spectrum utilized for the analysis. 

We normally combine the radial velocity measurements 
from spectra taken within a short time frame to give one 
velocity value. For VB 10, we binned all observations 
obtained over 2.5 hr intervals (4 - 16 individual measure- 
ments) to give the final data set used in this study be- 
cause this is a short amount of time relative to the orbital 
period of the putative planet we aimed to detect. The 
resulting radial velocities are given in Table [T] (with the 
best fit offset subtracted - see the next section) and plot- 
ted in Figure [TJ The typical uncertainty in the values is 
10m s -1 . 

3. ANALYSIS 
3.1. Comparison of models 

The radial velocity measurements we have obtained 
exhibit an rms dispersion of 10.6 ms -1 around a con- 
stant value when binned over 2.5 h, which is consistent 
with our estimated uncertainties. Therefore, we have 
not detected reflex orbital motion induced on VB 10 by 
an unseen companion. Therefore, the question becomes 
whether or not our data are inconsistent with the predic- 
tion of the existence of a giant planet from PS09, and the 
potential detection of radial velocity variability by Z09. 

Our measured radial velocities are plotted along with 
some possible orbit models in Figure [2] The data span 
82% of the orbital period (0.744 yr) suggested by PS09, 
in steps of ~ 0.25 phase units. Given that our uncertain- 
ties are roughly two orders of magnitude smaller than 
the semi-amplitude expected for a giant planet with an 
orbital period of a few hundred days (> lkms -1 ), it is 
likely that we would have detected significant changes 
in the star's radial velocity over the time span the data 
were obtained if such a planet exists. 

As discussed in §2, our measurements, like all other 
high-precision radial velocity measurements, are relative 
to an arbitrary scale. We have determined a best fit 
constant offset for all the orbits shown in this paper to 
account for this. We have fixed the orbital period and 
inclination (i — 96.9°) to the values given by PS09 for all 
the analyses presented in this paper because our aim is 
to assess their specific claim of a planet detection. The 
corresponding uncertainties in the period and inclination 
given by PS09 do not have a significant influence on the 
models relative to the precision of our data. For example, 
adjusting the assumed period and inclination within their 
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Fig. 2. — Top Measured radial velocities for VB 10 phased to 
the orbital period given by PS09 (points). The data are binned 
over each observing run for clarity and the error bars are smaller 
than the symbols. The dashed orange line gives the orbit model 
suggested by Z09 with the relative offset determined to give the 
best fit. The solid dark blue line gives the best fit to the data 
for a model with a 6Mj up planet on an orbit with e = 0.5. The 
light blue solid lines give the best fit models for 3Mj up (smaller 
semi-amplitude) and 9 Mj up (larger semi-amplitude) planets also 
with e=0.5. Bottom Residuals from the best fit 6Mj up planet 
model (points). The grey bar gives the range of values that would 
be consistent with three times the observed dispersion in the data 
(±32ms- 1 ). 



1 a uncertainties to obtain the smallest possible velocity 
semi-amplitude can only yield a maximum reduction in 
the semi-amplitude of < 3%. In addition to the period 
and inclination, we also assumed the same mass for VB 10 
that PS09 adopted (M* = 0.078 M Q ) to stay consistent 
with their study. 

The orbit model that was found by Z09 to give the 
best fit to both their NIRSPEC radial velocities and the 
PS09 astrometry is shown in Figure [2] for comparison 
with our data. The rms of the residuals of our data from 
this model with the relative offset optimized is 870 m s -1 . 
Therefore, the Z09 model is inconsistent with our mea- 
surements. 

Some additional orbit models for planets with masses 
of 3, 6, and 9Mj up are shown in Figure [2] This range of 
masses roughly corresponds to the nominal value and up- 
per and lower 1 a confidence limits determined by PS09. 
The time and location of periastron, and the eccentricity 
of the candidate planet detected by PS09 are not con- 
strained by their astrometry data. For these example 
models we adopted a fixed eccentricity e = 0.5, and op- 
timized the remaining parameters to give the best fit to 
the data (time of periastron, longitude of periastron, and 
radial velocity offset). As can be seen in the figure, the 
optimization of the unconstrained parameters leads to 
the orbital phases when the largest velocity variations 
are expected to be outside the times covered by our ob- 



servations. For reference, the best fit time of periastron 
is 2455103, and the best fit longitude of periastron is 
10°. However, the models do not fit the data even with 
these optimizations. The lower panel in Figure [2] shows 
the residuals from the M p = 6 Mj up model. The rms of 
the residuals is 256 ms , which is much larger than the 
observed dispersion in the data and the uncertainties in 
the individual measurements. 

3.2. Limits to a giant planet 

We conclude from the comparisons described above 
that it is unlikely that our radial velocity measurements 
are compatible with the presence of a giant planet around 
VB 10 with an orbital period of 0.744 yr and nearly edge- 
on orientation. Therefore, we were motivated to under- 
take a more specific and robust determination of what 
the data rule out. We did this by mapping the quality 
of orbit model fits to the data as a function of the pos- 
sible parameters, and comparing the determined values 
to the fit quality of a model representing no radial ve- 
locity variations (i.e. a flat line). It can be seen from 
inspection of Figure [5] that higher eccentricity and lower 
mass planets could be more easily hidden in the data. 
This led us to map the fit quality as a function of possi- 
ble planet mass and orbital eccentricity. We considered 
planet masses ranging from 0.5 to 15Mj up , and eccen- 
tricity values ranging from 0.0 to 0.95. We again fixed 
the period and inclination to those values given by PS09, 
and we marginalized over the unconstrained parameters 
at each grid point. 

Figure [3] shows the obtained map of the fit \ 2 - For ref- 
erence, the x 2 f° r a fit of a flat line to the data is 23.1 with 
11 degrees of freedom. This somewhat higher than ex- 
pected x 2 suggests that either VB 10 is exhibiting real ra- 
dial velocity variability at the 10ms -1 level, jitter arising 
from activity is affecting the measurements also at about 
the 10 ms -1 level (note that VB 10 is known to be a var i- 
able source with some flare activity, iBerger et al.| [2008'). 
we have underestimated our uncertainties by about 50%, 
or some combination of these effects. At this point it 
is not critical to understand the origin of the additional 
variability in the radial velocity measurements because 
this investigation is focused on ruling out much larger 
variations. 

The 5 a limit shown in Figure [3] corresponds to the 
boundary where the fit x 2 increases by more than 25 
from the x 2 of the flat line fit to the data. Everything to 
the lower right of this line is ruled out at high confidence. 
We find that for all conceivable eccentricity values (i.e. 
e < 0.95) we can rule out planets with M p > 3 Mj up and 
the orbital period and inclination suggested by PS09. 
This means that we can definitively rule out planets with 
masses in PS09's given 1 a confidence interval (indicated 
by the dotted lines in Figure [3]) . All higher mass planets 
are also ruled out. At the lower range of possible planet 
masses the determined limit is dependent on the assumed 
eccentricity. For example, planet masses of 1 Mj up and 
0.5Mj up can only be ruled out for e<0.7 and e<0.5 
respectively. 

4. DISCUSSION 

We have presented radial velocities for VB 10 that were 
obtained in the NIR using a new type of gas cell and 
that have a typical internal precision of 10 ms -1 . These 
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Fig. 3. — Map of the best fit \ 2 for a planet with P = 0.744 yr 
and i = 96.9°. At each point in the map the time of periastron, 
longitude of periastron, and radial velocity offset were optimized 
to give the best fit. The 5 a upper limit is given by the dashed 
contour. The region to the lower right of this line is excluded by 
our data. The dotted lines delineate the 1 a confidence interval on 
the planet mass from PS09. 



data were obtained in four epochs over a time span of 
0.61 yr. Our radial velocity measurements do not ex- 
hibit significant variability and are essentially constant 
at a level consistent with the measurement uncertain- 
ties. Therefore, we see no evidence for a giant planet 
orbiting VB 10 similar to that detected via apparent as- 
trometric perturbations by PS09. In addition, we do not 
confirm the ~ 1 km s _1 variability tentatively detected by 
Z09 using NIRSPEC radial velocity measurements that 
have a typical internal precision of 300 ms -1 . Our data 
are also inconsistent with the orbit model presented by 
Z09 as giving the best fit to both the PS09 astrometry 
and their radial velocities combined. If we consider our 
radial velocity data alone, then we can rule out plan- 
ets with Mp>3Mj up in the orbit suggested by PS09. 
Planets with masses down to 1 Mj up would need to have 
unusually large orbital eccentricities and be phased in 
a particular way to not have caused highly significant 
variations in our measured radial velocities. 

After we submitted our paper we became aware 
of another paper discussing the proposed VB 10 
lanet that was submitted at roughly the same time 
Anglada-Escude et al.ll2010l ). The authors of that pa- 
per presented radial velocity measurements for VB10 
with a typical precision of ^200ms _1 and that cov- 
ered roughly the same time period as our measurements. 
Their data also indicate a non-detection of the reflex mo- 
tion of VB 10 expected du e to the presence of th e giant 
planet proposed by PS09. lAnglada-Escude et all t)2010f ) 
carried out a combined analysis of the PS09 astrometry 
and their radial velocities and found that the possible 
orbit models do not give a significantly better fit to the 
data than a model assuming no reflex motion. They 
concluded that the giant planet proposed by PS09 is not 
likely to exist on the basis of these results, although they 
couldn't completely rule it out with their data. 

Our results seem to completely rule out the presence of 
a giant planet similar to that pro posed by PS09, which 



(|2010f ). Therefore, we suggest that the planet detection 
claimed by PS09 is spurious. If further work verifies this 
result then it would still be the case that no planet has 
yet been discovered via astrometry, and that the cur- 
rent state-of-the-art for the astrometric study of exo- 
planets remains post priori detection and characteriza- 
tion of a few select planets that we re first identified with 
the radial velocity method (e.g. [Benedict et al.l [20021: 
McArthur et all 12004 iBean et all l2007t iMartioli et al.l 



20101) . We note that obtaining astrometric measurements 



at the level of precision needed to detect exoplanets is 
notoriously difficult, especially from the ground as PS09 
have attempted to do. There is a long history of claimed 
and later refuted astrometric cxoplanct detections dating 
back to the famous case of the possible planet around 
Barnard's star, which was ori ginally suggest e d bas ed on 
25 years of obs ervations by Ivan de Kanrpl (119631) but 
later r efuted by iHershevl (|1973f ). iGatewood fc Eichhornl 
(|1973f) . and others. Furthermore, PS09 were not even 
the first to propose the existence of a sub-stellar com- 
panion t o VB 10 based on app arent astrometric pertur- 
bations (|Harrington et al.lll983|) . The first candidate was 
later disproved by subsequent astrometric m easurements 
(|Monet et alj|1992t lHarrington et al.lll993l ). 

With regards to the possible radial velocity variations 
in VB 10 seen by Z09, we suggest their analysis could 
have been influenced by an important systematic effect. 
This is a possibility that those authors themselves ac- 
knowledged. Of the five measurements presented by Z09, 
the one point in their data set that is inconsistent with a 
constant radial velocity (i.e. no variations possibly due 
to orbital motion) is also based on a spectrum that was 
taken with a different slit width than all the other data. 
Z09 did not use a deconvolution technique to account for 
changes in the spectrograph point spread function, so the 
velocity measured from this spectrum could have a differ- 
ent relative zero point compared to the other data. It is 
not clear whether the different slit width would cause as 
large a discrepancy as they saw (~lkms _1 ), but it is a 
possible candidate for reconciling the difference between 
their measurements and ours. 

Despite the negative result indicated by our data, we 
propose that there is a positive general result from this 
work. Namely, that we have demonstrated the power 
of NIR radial velocities obtained using CRIRES at the 
VLT with our ammonia cell. Late-type M dwarfs are 
an interesting potential planet host sample, but they are 
also extremely challenging for the normal planet search 
methods due to their faintness. For example, VB 10 is 
one of the closest, and therefore brightest, stars with an 
estimated mass < 0.1 Mq, but it still has V = 17.5. This 
puts it far beyond the reach of high-precision radial veloc- 
ity measurements in the visible with current instruments. 
However, late M dwarfs are reasonably bright in the NIR 
due to their extreme redness, and VB 10 has K = 8.8. By 
developing a technique for measuring high-precision NIR 
radial velocities, we have opened the door to a sensitive 
search for planets around stars at the very end of the 
main sequence and even beyond. The method is a valu- 
able stand-alone tool, and also a complement to the ongo- 
ing astrometric planet searches targeting M dwarfs (e.g . 
Pravdo et al1l200l iHenry et al.l [2001 iBoss et all [20091: 
Forhrich fc Bergen I2009D . We are continuing to obtain 
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is consistent with the results of lAnglada-E scude et"aT1 radial velocity measurements for VB 10, and many other 
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very low-mass stars to search for planets. The results 
from this work will help constrain the nature of planets 
around the most numerous stars in the Galaxy. 
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